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Abstract. We obtain accurate resonance energies for the Schrodinger equation with 
a central-field potential by means of a method based on a rational approximation to 
the logarithmic derivative of the wavefunction. We discuss the rate of convergence 
of our approach and compare present results with those obtained earlier by other 
authors. We show that present method is superior to the spherical-box approach 
applied recently to the same problem. As far as we know present results arc more 
accurate than those available in the literature and may be a suitable benchmark for 
testing future approaches. 



1. Introduction 

In a recent paper Zhou et al [T] applied the well-known spherical-box stabilization 
method to the calculation of the resonance energies of the Schrodinger equation with 
the potential V(r) = Vor 2 e~ r + Z/r. They integrated the eigenvalue equation by means 
of the Runge-Kutta method and estimated the positions and widths of the resonances 
from the behaviour of the bound-state energies as functions of the box radius. In 
particular they calculated the first two s-wave the first p-wave and the first d-wave 
resonances. Zhou et al [1] experienced some difficulties in estimating the position an 
width of the second s-wave resonance and could not obtain the third one. They also 
obtained rather crude estimations of the first p-wave and d-wave resonances. It seems 
that the spherical-box approach is rather ill-suited to broad resonances. 

The potential V(r) mentioned above has proved a suitable benchmark for the 
development and testing of several methods for the calculation of the energies of 
metastable states pH2l]. Most authors have considered the case Z = [2TT21] and 
just a few ones included the Coulomb interaction Z — — 1 [20} 123]. 

Some time ago, we applied the Riccati-Pade method (RPM) to the calculation of 
the lowest s-wave resonance of the potential V(r) with Z = [17]. In that earlier paper 
we did not discuss the rate of convergence of the method on this particular model and 
merely showed the result for the lowest s-wave resonance. The purpose of this paper 
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is to compare the RPM with the spherical-box approach and with other alternative 
methods on the central-field potential with Z — — 1. 

In Sec. [2] we outline the method. In Sec. [3] we apply the RPM to the potential V(r) 
with Z = — 1, analyze its results and compare them with those obtained earlier by other 
authors. Finally, in Sec. H]we discuss the advantages of the RPM and draw conclusions. 



2. The method 

The radial part of the dimensionless Schrodinger equation for a central-field potential 
V(r) is 

1 d 2 1(1 + 1) 



V(r) 



$( r ) = E${r) (1) 



2 dr 2 2r 2 

where I = 0, 1, ... is the angular-momentum quantum number and 3>(0) = 0. The RPM 
is based on a rational approximation to the regularized logarithmic derivative of the 
wavef unction 

m = i±i _ m (2) 

r <P(r) 
that can be expanded as follows: 

oo 

f(r) = E & (3) 

3=0 

Note that the term (/ + l)/r removes the singularity of $'(r)/$(r) at origin and 
that we can obtain the coefficients fj(E) analytically by means of simple recurrence 
relations [T7] . 

We then convert the Taylor series into a rational approximation: 

s T N+d n ri 2N+d+i 
i N + d/m = ^ L 1 f-= E ff 3 (4) 

where N — 1,2, . . . and d — 0,1, . . .. The 2N + d+ 1 adjustable coefficients aj and bj are 
insufficient to provide the 2^ + ^ + 2 coefficients fj. This condition is satisfied only if the 
Hankel determinant Hf)(E) with matrix elements fi + j + d-i{E), i,j = l,2,...,D = N+l, 
vanishes [17] (and references therein). The RPM conjecture is that there are sequences 
of roots EPA of H d D (E) = 0, D = 2, 3, . . . that converge towards the actual bound- and 
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metastable-state energies of the Schrodinger equation (pQ). The calculation is remarkably 
simple because the Hankel determinants are polynomial functions of the energy. 



3. Results and discussion 

For comparison we consider the potential 

V(r) = V r 2 e~ r + - (5) 

with the model parameters Vo = 7.5 and Z = —1 [IJ[20j[23]. We label the complex 

energies of the metastable states Ei >v so that Re E^+i > ReEi yU , v — 0, 1, The s-, 

p- and d-waves discussed by Zhou et al [1] and Sofianos and Rakityansky [23J correspond 
to I — 0, I — 1, and / = 2, respectively. 

Since we are looking for Siegert pseudo states that satisfy [25] 

lim —±4- = ik (6) 

then it seems reasonable to choose d = because 

«m[W/JV] = £ (7) 

For that reason it should be assumed that d = from now on, unless otherwise stated. 

In order to estimate the rate of convergence of the RPM we calculate = 
log |a£> — «_d + i| where old is either the real or imaginary part of E^ D '°K Fig. [1] shows that 
the rate of convergence of the RPM for both the positions and widths of the first three 
s-wave resonances is remarkable. It is worth mentioning that in the case of a narrow 
resonance the imaginary part of the root will appear at sufficiently large determinant 
dimensions D; that is to say, when | HeE^ D ' ^ — Re E^ D+l ^ \ is of the order of magnitude 
of | Im E\ . We appreciate that the rate of convergence (given approximately by the 
slope of Ld) is almost independent of u; the main difference is that the greater the 
value of v the larger the determinant dimension D necessary for the appearance of 
the corresponding sequence. On the other hand, the performance of the spherical-box 
approach deteriorates as the resonance width increases [I]. It is clear that the RPM is 
preferable to the spherical-box approach, at least for this example. 
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Fig. [2] shows that the rate of convergence of the RPM for the p-wave resonances 
(I = 1) is similar to that discussed above. It is clear that the rate of convergence of the 
Hankel sequences is also independent of /. We confirm our conclusion that the RPM is 
preferable to the spherical-box approach because Zhou et al pQ roughly estimated the 
position and width of -E^o and were unable to obtain other p-wave resonances. 

Fig. [3] shows the rate of convergence for the first two d-wave resonances. The 
behaviour is similar to those discussed above for the s and p ones. According to Zhou 
et al [I] the spherical-box approach only revealed the first d resonance for which they 
could provide a rather crude estimate of the position and width. 

As far as we know, the most accurate results for this model are those calculated 
some time ago by Sofianos and Rakityansky [23] . Present results are even more accurate 
and may therefore be a useful benchmark for other approaches. For that purpose we 
show them in Table [D 

Finally, we mention that the rate of convergence of the RPM is not affected by 
the choice of the displacement d. In the present case, for example, we obtained similar 
results with d = 1 that we do not show here. 

4. Conclusions 

One of the main advantages of the RPM is its remarkable simplicity. We first obtain the 
coefficients of the Taylor series ([3]) by means of a straightforward recurrence relation |17j . 
Second, we construct the Hankel determinant, which is a polynomial function of the 
energy, and find its roots. Third, we identify the sequences of roots that converge 
towards physically acceptable results. 

Another advantage of the RPM is that exactly the same Hankel determinant applies 
to both the bound states and resonances. It comes from the fact that the RPM does not 
take explicitly into account the asymptotic form of the wavefunction at infinity and the 
rational approximation applies to any solution of the Schrodinger equation. Of course, 
we have to take into consideration the behaviour of the wavefunction at origin in order 
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to remove any singularity of $'(r)/$(r). 

We think that present results clearly show that the RPM is much more accurate 

and reliable than the spherical-box approximation. We have calculated the resonances 

discussed by Zhou et al p] with much more accuracy and also obtained others that those 

authors were unable to locate. Besides, it is worth noting that the RPM is as simple, 

or even simpler, than the box-stabilization method in any of its forms p^[6tlT3|[T¥j. 
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Table 1. Some complex energies Ei M for the potential ([5]) 
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1.7805245363623048 


0.00004785969842876 
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4.101494946209 


0.578627213766 
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4.6634610967 


2.6832007703 
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3.848001634811759 


0.137692229585768 
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4.750053489274 


1.75278992436148 
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4.9005161468291143 


0.7837535082665858 
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1 


5.3006134902578 


2.942357430621 




Figure 1. Convergence rate Lr> for the real (solid line) and imaginary (dashed line) 
parts of the energies Eq^ 
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Figure 2. Convergence rate Ld for the real (solid line) and imaginary (dashed line) 
parts of the energies E\ tV 




Figure 3. Convergence rate Lr> for the real (solid line) and imaginary (dashed line) 
parts of the energies i?2,i/ 



